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New boundary conditions are imposed on the familiar cylindrical gravitational wave vacuum 
spacetimes. The new spacetime family represents cylindrical waves in a flat expanding (Kasner) 
universe. Space sections are flat and nonconical where the waves have not reached and wave 
amplitudes fall off more rapidly than they do in Einstein-Rosen solutions, permitting a more regular 
null inifmity. 
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This paper shows how to construct exact solutions for compact cylindrical gravitational wave pulses spreading in an 
expanding universe that is spatially flat at large cylindrical radius, where the waves have not reached. The universal 
expansion is restricted to the translation-direction (or z-axis) of cylindrical symmetry. The remaining directions 
show gravitational waves spreading in a three-dimensional spacetime that can be exactly Minkowskian before the 
^vq , wave pulse. The effective three dimensional spacetime has a completely regular conformal compactification so that 
a radiation zone at future null infinity is well defined without any rescaling of fields. This simple behavior makes 
fS | these solutions ideal for testing the numerical simulation codes that are used to predict gravitational waves from 
compact sources. Because they have the simplicity of Einstein-Rosen solutions but are compatible with the boundary 
conditions of the real universe, these new solutions may model localized inhomogeneities in Big Bang cosmologies and 
may also provide an interesting restricted framework for testing formal ideas such as quantum gravity. 
' The original Einstein- Rosen cylindrical wave solutions of Einstein's vacuum field equations [l|, [2| inhabit a universe 
that is conical rather than flat at infinite cylindrical radiusQ and which, apart from the oscillations of the waves, 
■ does not evolve in time. Thus, they are not good models for inhomogeneities in our own universe, which is thought 
\& , to be spatially flat and expanding and do not make convenient test cases for numerical simulations of gravitational 
i waves from compact sources. However, the boundary conditions can be changed. Cylindrically symmetric (G2) exact 
solutions with closed space sections of S 3 , S 1 x S 2 and S 1 x S 1 x S 1 (three-torus) topologies are known and have 
proven useful as toy models of inhomogeneous cosmologies. [3, [f| 0] The three-torus solutions, in particular can be 
regarded as standing-gravitational- wave inhomogeneous generalizations of spatially flat Kasner universes 0J1> S E3] 
and have provided useful test cases for validating numerical simulations of Einstein's field equations. Ill, HI] Here, 
^ | we will change the boundary conditions again. 

il We will use the notation of earlier papers[H [U, Q on these spacetimes, but with some changes in coordinate names 
to fit the new circumstances. The isometry group coordinates are z and ip (instead of a and S used in the earlier 
rr| papers), with z ranging from — 00 to +00 and cp ranging from to 2tt. The remaining coordinates t,r each range 
^ ■ from to 00. The isometry group consists of translations of z and of tp modulo 2ir. Thus, each group orbit is an 
intrinsically flat cylinder. 

The advanced and retarded time coordinates, u = t + r and v = t — r, govern the propagation of waves in these 
geometries. It is convenient to use subscripts to denote derivatives with respect to advanced and retarded time so 
that, for example, i?+ = dR/du and i?_ = dR/dv. The spacetime metric in these coordinates is 

ds 2 = e 2a (dr 2 - dt 2 ) + R (e 2 ^dz 2 + >d<p 2 ) . (1) 

Two of the functions in this form of the metric satisfy wave equations, 

R+- = 0, (2) 

and 

+ (Jty_)+ = 0, (3) 
while the remaining function, a, is determined by integrating the equations 

R+a + - R^l + ±R + - ±5± (4) 



and 



i?_a_ = Rip 2 _ + - -J^-. (5) 
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Because the area of a group orbit is proportional to R, the global character of the solution is established by choosing 
a particular solution for R. The Einstein-Rosen cylindrical wave solutions correspond to the choice R = r. For the 
solutions sometimes referred to as Gowdy T3 models, the choice is R = t. For spacetimes with space sections of S 3 
or S 1 x S 2 topology, the appropriate choice is R — sinr sin 0,0, Q Here, we explore a new choice: 

R = rt. (6) 

All of the orbits have zero area at t — and also at r = 0. This spacetime expands from an initial singularity at t = 
and has a cylindrical symmetry axis at r — 0. 

Separate a logarithmic solution of the wave equation from the function i/> 

ip = ~hi(r/t) + W (7) 

so that the metric becomes 

ds 2 = -e 2a dt 2 + e 2a dr 2 + e~ 2W r 2 d V 2 + e 2W t 2 dz 2 , (8) 
where W satisfies the same wave equation as ip. There will be a conical singularity unless the function a obeys 



da 
dr 



0, (9) 



and the ratio of coefficients of the radius term dr 2 and the angle term r 2 d(p 2 is equal to its flat space value of 1 to 
second order in r. The result is a non-conicality constraint, 

(a + W% =0 = 0, (10) 

and a derivative constraint on the function W. 



dW 
dr 



0. (11) 



The following general argument, which depends on the presence of an initial singularity, shows that there should 
exist solutions that are spatially flat at large cylindrical radius: Because W satisfies a wave equation, its value on 
a spacelike hypersurface E near the initial singularity may be specified arbitrarily. In particular, W and its time 
derivatives can be specified to have compact support S on the hypersurface so that they vanish everywhere else on 
£ and will therefore vanish throughout the corresponding future domain of dependence D + (E — S). 13] Wherever 
W is zero, Einstein's equations imply that the function a satisfies a homogeneous wave equation and may then be 
zero throughout this same domain of dependence. Thus, it should be possible to have a compact wave source in a 
spacetime whose metric at large distance from the source has the Flat Kasner formflij 

ds 2 = -dt 2 + dr 2 +r 2 dip 2 + t 2 dz 2 , (12) 

which is the result of setting a = W = in ([8]) . (This form of the metric is discussed further near the end of this 
paper.) The constant-z surfaces in such a spacetime are intrinsically flat and do not show the conical behavior that 
is familiar from the Einstein- Rosen cylindrical wave solutions. In order to achieve this regularity, it must be possible 
for the function a to satisfy a second non-conicality constraint; aL^^ = 0. 

Now construct such a solution. From © and 0, the function W obeys the wave equation 

1 d dW 1 d r dW _ Q 
t dt dt r dr dr 



for a scalar field in the Flat Kasner metric (|12[) above, with solutions 

W = { dk[A (k) J (kt) + C (jfe) iVo (Art)] J a (kr) , (14) 



where Jo and A^o are Bessel functions of the first and second kind. The remaining metric function, a, is determined 
by integrating equations J3J and (J5j) , which now become 

ua+ = - (u 2 - v 2 ) W\ - vW+ (15) 
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and 

-va- = i (u 2 -v 2 )W 2 +uW^. (16) 

This last equation implies a matching constraint at v — or r = t of the same sort that is encountered in the S 3 and 
S 1 x S 2 solutions. [J, 

rW-\ r=t = Oi-1. (17) 

Note that the quadratic nature of (|16|) yields two solutions and thus two possible forms of the matching constraint. 
This constraint is preserved in time when the field equations are satisfied. Impose it on the solution given by (|14J) and 
obtain an expression that involves only the Wronskian of the Bessel functions. Abel's formula for the Wronskian jl5l [lr| 
then yields a constraint on the coefficients C (k) : 

dk 1^1 = o ;7T . (18) 
k 

Here, we can see that one possible form of the matching constraint permits C (fc) = while the other requires 
C (k) ^ 0. When either form of this constraint is satisfied, the function a can be expressed as the integral 

f x t 2 Z 

a(x,t) + W(x,t) = I t2 _ r2 rdr, (19) 



where 



Z = (W 2 + W 2 ) - 2-W r W t - 2^W r + -W t . (20) 

In the case of interest, with W (x, t) = for x 3> t, the requirement of a non-conical metric in the part of spacetime 
not reached by gravitational waves means that there is an integral constraint on the initial data for the wave: 

oo ,27 

- 2 - 2 rdr = 0. (21) 

q t z — r z 

This constraint is quite similar in form to the integral constraints that are encountered in the spatially closed G2 
solutions. In those cases, the constraint is needed to have a metric that is non-conical at both cylindrical symmetry 
axes or, for the three-torus solutions,a metric that is periodic in space. Here, one of the 'axes' is at infinity but the 
situation is otherwise the same. 0, Q 

Notice that for t > r the function Z is dominated by the positive definite expression (W 2 + W 2 ) and can be 
thought of as an energy density. In that region, the solution is just a standard Einstein-Rosen cylindrical wave in 
distorted coordinates. At large values of t, the coordinates go over to the canonical Einstein-Rosen coordinates and the 



integrand becomes what has sometimes been called the "C-energy" . [17} The positive definite nature of the C-energy 
is the reason that standard Einstein-Rosen cylindrical wave solutions must be conical at large cylindrical radius. 

In the region where r > t, the integrand reverses sign and the dominant term of Z becomes W r Wt, which has the 
form of a scalar field energy-flux and can have either sign. In this region, our solution is just a distorted form of the 
cosmological Einstein-Rosen solution (sometimes identified as a Gowdy T 2 or T 3 toroidal spacetime) that expands 
from an initial singularity. In the canonical coordinates for this type of solution, the analog of our integral for the 
function a contains only an energy-flux type of term. The form of the integral constraint makes it clear that there 
must always be some non-zero wave in the r > t region. A wave with compact initial data near r = t = would 
straddle the r = t surface with parts in both regions. A wave pulse with compact initial data near r = at a later 
time would need a precursor pulse with r > t to satisfy the constraint. 

A well known disadvantage of Einstein-Rosen cylindrical waves is that their amplitudes fall off as r -1 / 2 at large 
cylindrical radius r, which prevents them from being regular at future null infinity. The difficulty can be managed 
by defining rescaled wave amplitudes. [l8j but that special treatment complicates the use of cylindrical waves as test 
cases for formalisms and numerical methods designed to treat three dimensional waves. The solutions discussed here 
acquire an additional factor of t~ x / 2 from their Bessel function time dependence so that wave packets with support 
near the r = t surface fall off as 1/r. 

A simple way to examine future null infinity is to apply the compactified coordinates, r and 9 defined in the 
usual way by u — tan i (r + 8) and v — tan i (r — 9) to the three dimensional spacetime perpendicular to the 



84014 



symmetry axis. [HI l20j The metric 3 ds 2 of this spacetime is brought into the form 3 ds 2 — Q 2 ds 2 where f2 = 
2 cos | (r + 6) cos \ (r — 9) and 

dS 2 = e 2a (d6 2 - dr 2 ) + sin 2 6e- 2W d^\ (22) 

so that the spacetime is embedded as the open set defined by0<r<7r, 0<r + #<7rma cylindrical S 2 x S 1 
spacetime. Future null infinity is then the set of points described by r + 9 = tt. For Einstein-Rosen waves one finds 
that W has the form W — ► fJ 1 / 2 / where / is a smooth function, so that derivatives of W will contain Or 1 / 2 and will 
diverge at future null infinity. Thus, Einstein-Rosen waves do not admit a regular compactification of the spacetime 
perpendicular to the symmetry axis. For the solutions discussed here, W has the form W — > f2 / and is clearly regular 
at future null infinity. Because the derivatives of the function a are related to the derivatives of W by (TT5]) and (flB|) 
that function is also regular at future null infinity. The solutions presented here have a regular compactification with 
no need for rescaling. 

A more complete analysis of null infinity in Einstein-Rosen spacetimes indicates that the full four-dimensional 
spacetime is somewhat better behaved than the three-dimensional section perpendicular to the symmetry axis.[2ll.l22l| 
It would be interesting to see a similar analysis of the solutions presented here. It is worth noting that null infinity 
for the toroidal model spacetimes has also been analyzed. [23l l2~ij 

The Flat Kasner metric given by (|12[) serves as the background spacetime for these gravitational waves. Here, we 
choose to regard it as an axially symmetric Bianchi type I homogeneous spacetime expanding from a singularity at 
t = 0. However, it can be characterized in several other ways that may be significant for a complete analysis of null 
infinity in these solutions. It has zero spacetime curvature and is similar to the Milne cosmology. [25| It is the sector of 
Minkowski spactime that corresponds to the future of the t = z = two-surface and is presented in Rindler hyperbolic 
coordinates. The symmetry in our z direction corresponds to Lorentz boosts along one Minkowski space axis. If our 
z coordinate is made to be periodic, the result is Misner space[26]. With additional periodicities, the result is the 
background geometry of the toroidal Gowdy T3 spacetimes. The fact that the Milne and Misner spaces have long been 
favorite toy models of Big Bang universes adds to the possible usefulness of our solutions. 

The solution for the wave function W given by (|14[) may not always be the most convenient format to use. For 
example, one might wish to give Cauchy data at an initial time instead of amplitude functions. That problem has 
been solved by Woodhouse for the conventional Einstein-Rosen waves and a similar analysis may be possible in 
this, somewhat different case.[27[ Again, one might wish for a less general set of solutions that are given in terms of 
elementary functions so that one can avoid the need to perform integrals. That problem was solved for the conventional 
Einstein-Rosen waves by the Weber- Wheeler pulse solution and a similar specialization should be possible here as 
well. |] 

The solutions presented here have just one polarization of gravitational waves and no matter fields. It is well known 
that a second polarization can be added in the form of an off-diagonal term in the metric Q and massless matter fields 
such as an electromagnetic wave [28| and a massless scalar or axion field (2^, [3(| can also be added without changing the 
basic structure of these spacetimes. The effect of the added degrees of freedom is to replace the linear wave equation 
([5]) or (|13|) by a system of nonlinearly coupled wave equations. These nonlinear generalizations have provided a simple 
setting for numerical explorations of nonlinear gravitational interactions. [{j, [qL floT| Analytic procedures for genera ting 
solutions to this nonlinear system, as it appears in canonical Einstein- Rosen coordinates, have also been developed. [3l| 

Because the Einstein- Rosen solutions and their nonlinear generalizations are the simplest examples of gravitational 
waves that radiate outwards, they have been used to test both formal ideas and numerical procedures in a simple 
set ting . There is, for example, a long history of using them as a framework for testing ideas about quantum gravity. [30l 
I321I33I] Thev have, as another example, been used to develop Cauchy characteristic matching techniques for extracting 
gravitational radiation from numerical simulations of astrophysical sources. [TH |34| Another technique to accomplish 
the same purpose, the use of hyperboloidal slicing, is still being tested on scalar waves in fixed background spacetimes 
but will need more complex tests shortly. [35| The solutions presented here are slightly more complex than the Einstein 
Rosen solutions. However it may be argued that they fit the boundary conditions of our own universe while the 
Einstein Rosen solutions do not. Thus, these solutions, and their nonlinear generalizations to multiple polarizations 
and massless fields are a logical step up from the Einstein Rosen solutions as test vehicles for new techniques. 

Note added in proof. - The solutions presented here belong to the general class of boost-symmetric spacetimes 
discussed by Bicak and Schmidt, Phys. Rev. D 40, 1827 (1989). 
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